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Abstract. For the last one and a half decades it has been known that the 
exponential product formula holds also in norm in nontrivial cases. In this 
note, we review the results on its convergence in norm as well as pointwise of 
the integral kernels in the case for Schrodinger operators, with error bounds. 
Optimality of the error bounds is elaborated. 
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1. Introduction 

The Trotter product formula, Trotter-Kato product formula or exponential prod- 
uct formula is usually a product formula which in strong operator topology ap- 
proximates the group/semigroup with generator being a sum of two operators. It is 
often a useful tool to study Schrodinger evolution groups/semigroups in quantum 
mechanics and to study Gibbs semigroups in statistical mechanics. 

To think of a typical case, let A and B be selfadjoint operators in a Hilbert 
space TC with domains D[A] and D[B] and H A + B their operator sum with 
domain D[H] = D[A]r]D[B]. Assume that H is selfadjoint or essentially selfadjoint 
on D[H] and denote its closure by the same H. Then Trotter 03] proved the 
unitary product formula 



■^-itB/2n^-itA/n^-itB/2n^n 
I —itA/n —itB/nyn 




0, strongly, 

0, strongly, n 



oo. 



2 



Takashi Ichinose and Hideo Tamura 



and also, when A and B are nonnegative, the selfadjoint product formula 

[g-tB/2„g-tA/„g-tB/2„]„ _ ^-tH ^ 0^ strongly, 

[e-*A/n^-tB/n^n _ ^-tH ^ strongly. 71 ^ oo. 

The convergence is locally uniform, i.e. uniform on compact ^-intervals, respectively 
in the real line R and in the closed half line [0, oo). Kato [221 dicovered the latter 
selfadjoint product formula to hold also for the form sum H := A+B with form 
domain D[H^^^] = D[A^^^] n D[B^^^], which we assume for simplicity is dense 
in Ti. However, it remains to be an open problem whether the unitary product 
formula for the form sum holds. 

However, since around 1993 we have begun to know that selfadjoint product 
formulas converge even in (operator) norm, though in some special cases, by the 
following two first results. Rogava [37] proved, when B is A-bounded and H — 
A + B is selfadjoint, among others, the abstract product formula that 

||[e~*'^/"e-*^/"]"-e-*-f^|| = Oin"^^^ logn), n -> oo, 

locally uniformly in [0, oo). Helffer 13 proved, when H := —A + V{x) is a 
Schrodinger operator in L^(R'') with nonnnegative potential V{x) satisfying |9"F(a:)| 
< Ca {\a\ > 2) so that H is selfadjoint on the domain D[—A]r\D[V], the symmetric 
product formula that 

II |g-ty/2„g-t(-A)/ng-ti//2„]„ _ g-tff II ^ o{n-^), n ^ oo, 

locally uniformly in [0, oo). Many works were done to extend these results before 
2000, e.g. in [5llMl22l[3l[33l[35] for the abstract product formula, [i[T0l[17l[18l[19l 
il] for the Schrodinger operators, and after that, e.g. in [23l [27l [16] . [3ll4l[6l[7] for 
the abstract product formula. In most of them, use was made of operator-theoretic 
methods, though of a probabilistic method in [TT ] fTS ] [T9 l |4T]. 

In this note, we want to describe more recent results on convergence in norm 
for exponential product formulas and also pointwise of the corresponding integral 
kernels, mainly based on our works since around 2000, [231 [571 [Ml [HI [2S]- As for the 
error bounds, although it is easy to see by the Baker-Campbell-Hausdorff formula 
(e.g. |45j . [40j ) that with both operators A and B being bounded, the nonsymmetric 
product formula has an optimal error bound 0{n~^) while the symmetric one does 
0(n~^), it was shown in [23 that even the symmetric product formula has an 
optimal error bound 0{n~^) in general, if both A and B are unbounded. However, 
in [25] (cf. [26]), a better upper sharp error bound 0{n~^) has been obtained for 
the symmetric product formula with the Schrodinger operator — A + V{x) having 
nonnegative potentials V{x) growing polynomially at infinity, in spite that both 
—A and V are unbounded operators. In this note we mention, with a sketch of 
proof, a latest complementary result [2] which settles the sharp optimal error bound 
is in fact 0(n~^) with the symmetric product formula for the harmonic oscillator, 
by estimating the error not only from above but also from below, in norm as well 
as pointwise. 
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Theorems are described in Section 2. Optimal! ty of error bounds is discussed 
separately in Section 3. The idea of proof is briefly mentioned in Section 4. In Sec- 
tion 5 we give concluding remarks, and also refer to a connection of the exponential 
product formula with the Feynman path integral. 

It should be also noted that in almost the same context with the notion of 
norm ideals (e.g. [12], [38 ) we are able to deal with the trace norm convergence as 
in [461 EOl [311 [111 [14, 21, 42 . For an extensive literature on this we refer to [47 . 

The content of this note is an expanded version of the lecture entitled "On 
converegence pointwise of integral kernels and in norm for exponential product 
formulas" given by T.I. at the International Conference "Modern Analysis and 
Applications (MAA 2007)", Odessa, Ukraine, April 9-14, 2007, which is a slightly 
extended version of of the lecture (unpublished) given at the Conference on ''Heat 
Kernel in Mathematical Physics" , Blaubeuren, Germany, November 28-December 
2, 2006. 

2. Theorems 

We begin with our result which extends ultimately Rogava and Helffer's. 

Theorem 2.1. (Ichinose-Tamura-Tamura-Zagrebnov 2001 \'2'Sl \Tl\ ) Let A and B he 

nonnegative selfadjoint operators, and assume H = A+B is selfadjoint on D[H] = 
D[A] n D[B]. Then as n -> oo, 

||jg-ti3/2„g-tA/„g-ts/2„j„_g-tff|| ^ O(n-i), (2.1) 
||jg-tA/«g-tB/„]„_g-tff|| ^ 0(„-i). (2.2) 

The convergence is locally uniform in the closed half line [0,oo), while on the whole 
half line [0,oo), if H is strictly positive, i.e. H > rjl for some 77 > 0. The error 
bound 0(n~^) in h2.1]) and h2.1]) is optimal. 

We can go beyond this result. First, focussing on the Schrodinger operator 
—A + V{x), we ask whether norm convergence implies pointwise convergence of 
integral kernels. The answer is yes, though strong convergence does not. This prob- 
lem is discussed for Schrodinger operators with potentials of polynomial growth 
fTheorem 12. 2p . with positive Coulomb potential (Theorem 12.31) . and also for the 
Dirichlet Laplacian (Theorem 12. 4p . Pointwise convergence of integral kernels for 
Schrodinger semigroups is important, because it gives a time-sliced approximation 
to the imaginary-time Feynman path integral. 

Next, we ask, for the unitary exponential/Trotter product formula, whether 
there are nontrivial cases where it converges in norm, though it does not in general 
hold (see [H]). The answer is yes. In fact, it holds for the Dirac operator and 
relativistic Schrodinger operator (Theorem 12. 5p . 

Let H = Hq + V :— —A -I- V{x) with V{x) a real-valued function. By 
K'-"\t,x,y) we denote the integral kernel of [e-tHo/2n^~tv/n^-tHo/2n-^n^ ^^^^ 
e'*"{x,y) that of e~*^. 
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Theorem 2.2. (Ichinose-Tamura 2004 [25 j (positive potential of polynomial growth) 
Assume that V{x) is in C°°{'R.'^), bounded below and satisfies \df^V{x)\ < Ca{x)"^~^^ 
with some <S <1 ({x) = {1 + x^fl"^). 
(i) (In norm) 

||jg-tffo/2ng-ty/«g-tHo/2n]n _ ^-tH^^^^ ^ 0(?^"^), (2.3) 

locally uniformly in the open half line (0, c»). 
(a) (Integral kernel) 

[K^-\t,x,y)^e-'"{x,y)]^Oin-% 

in C°°(R'' X R'')-topology, locally uniformly in (0,oo), (2.4) 

i.e. together with all x, y-derivatives. 

This theorem improves the result of Takanobu [H] , who used a probabilistic 
method with the Feynman~Kac formula (see Sect. 5) to show uniform pointwise 
convergence of the integral kernels, roughly speaking, with error bound 0{n-p/^), 
if V{x) satisfies V{x) > C(l + |a;|2)^/2 and \d^V{x)\ < Ca(l + |x|2)(p-«I"I)+/2 for 
some constants C, Cq > and p > 0, < 6 < 1. The claim of Theorem 2.2 is a 
little bit sharpened in Theorems 3.1 and 3.2, in the next section, in the case of the 
harmonic oscillator. 

Theorem 2.3. (Ichinose-Tamura 2006 |26j ) (positive Coulomb potential) Let H := 
-A + V{x) with V{x) > 0. Assume that V{~A + 1)'": L^(R'') ^ L^(R'') is 
bounded for some < a < 1, and that V G C°° near a neighbourhood U of both p 
and q (after p,q G R'* taken). Then 

[Kit/nnx,y)^e-'"ix,y)]^Oin-'), 

in C°°(t/)-topology, locally uniformly in (0, oo). (2-5) 

The condition is satisfied if is in L"^ (R?) + L°° (R.^) , in particular, iiV is the 
positive Coulomb potential l/\x\. We don't know what happens at the singularities 
of y(a:). 

Theorem 2.4. (Ichinose-Tamura 2006 f26]j (Dirichlet Laplacian) Let fl C R'* be 

a bounded domain with smooth boundary and Xfi the indicator function offl. Let 
Ho := — A in L'^(Rf'), and H :— — A^ the Dirichlet Laplacian in with domain 
D[H] = H^{Q) n H^[n). Then /or < cr < i, 

(Xa e-'""/" xn^ix, y) - e-*^(x, y) = 0(n--), 

locally uniformly in {t,x,y) G (0,oo) x fl x fl. (2.6) 

We don't know what happens when x ot y approaches the boundary of fi. 

Corollary. 

II [xn e-*^«/" XnTf - e-*«/|U. ^0, / G L^n). 



Consequently, Theorem l2.4l is a stronger statement than this corollary, though 
the latter is also obtained by Kato as an abstract result: If A is a nonnegative 
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selfadjoint operator and P an orthogonal projection in a Hilbert space H, then 
(Pe-*^/"P)" -> e-*^^, strongly, as n ^ oo, where Ap {A'^/'^P)*{A^/'^P). In 
passing, however, it is an open question whether it holds that (Pe~**^/"-P)" 
^-itApp^ strongly {Zeno product formula). A partial answer was given in 

All Theorems 2.2-2.4 hold with order of products exchanged, e.g. in Theo- 
rem[221 [e-*^/2"e~*^°/"e-*^/2n]« instead of [e-*^o/2"e-*^/"e-*^«/2"]". 

Theorem 2.5. (Ichinose-Tamura 2004 [24\) (Unitary Trotter in norm) Let A and 
B he selfadjoint, and assume H :~ A + B to he essentially selfadjoint in a Hilhert 
space Ti.. Assume that there exists a dense suhspace V ofTi. with V C D[A] nD[B] 
such that e"'*-^, e~"^ : V V. Further assume that the commutators [A,B], 
[A, [A, B\\ and [B, [A, B]] are hounded on TL. Then 

||(e-"^/2"e-^*^/"e-^*^/2")"-e-"^|| =0(n-2)^ n oo, (2.7) 
locally uniformly in the real line R. 

As important applications we have ones to the Dirac operator H — Hq + 
V = {ia ■ V + m/3) + V{x) in ^^(R^)'', where a = (ai,a2,a3) and /3 are the 4 
Dirac matrices, with 9^y(a;)(|7| = 2) being bounded, as well as to the relativistic 
Schrodinger operator H = Ha + V = V-A + + V{x) on ^^(R'') with d2V{x) 
being bounded for 1 < I7I < 4 (0 < I7I < 4, if m = 0). In these cases, H are 
essentially selfadjoint, and satisfy the conditions in the theorem. So it holds that 

II ^^-^tV/2n^-^tHo/n^~^tV/2n^n _ ^-^tH^^^^ ^ ©(n'^), n ^ OO, (2.8) 

locally uniformly in R. 

However, this theorem does not apply to Schrodinger operators except for 
the Stark Hamiltonian (—A + V{xy) + a • x in Lp'iTL'^), where a is a constant real 
vector in R'^. 

Finally it should be noted that we have shown in Theorems 2.2-2.4 that the 
convergence is unifrom only on compact t-intervals which are away from 0, though 
in Theorems 12.11 and 12.51 on ones which are allowed to be not away from 0. 

3. Optimalty of Error Bounds 

In this section we discuss optimality of error bounds. The error bound 0(l/7i) in 
Theorem 12.11 is optimal, because if both A and B are bounded operators, by the 
Baker-Campbell-Hausdorff formula we know 

for some R'^ and i?„ being uniformly bounded operators which in general are not 
the zero operator. From this, optimality in the former non-symmetric case is evi- 
dent. But even in the symmetric case it is optimal. Indeed, there exist unbounded 
nonnegative selfadjoint operators A, B such that H = A + B is selfadjoint and 

||[g-*B/2„g-tA/„g-tS/2n]n _ g-*H|| > 
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for some continuous function c{t) with c{t) > 0, t > and c(0) — f[27j). 

However, further in some special symmetric case in Theorem l2 . 21 where — A, V 
are taken as A, B, we have seen the symmetric product formula hold with a sharp 
error bound 0{n^^). We can make more precise this result with the 1-dimensional 
harmonic oscillator H := Hq + V -.^ ^("^x + ™ i^(R). 

Theorem 3.1. (Azuma- Ichinose 2007 [2\) There exists bounded continuous func- 
tions C{t) > and c{t) > in t > 0, which are positive except t — with 
C(0) = c(0) = 0, independent of n, such that for n = 1,2,..., 

c{t)n-^ <\\[e-^^^e-^"°e-^^]" ~e~'"\\<C{t)n-^, t>0. (3.1) 

This theorem mentions an error bound from below, extending the harmonic 
oscillator case of Theorem 12.21 which treats only the right-half inequality with 
C(t) = C being a positive constant depending on each compact i- interval in the 
open half line (0, oo). 

It is anticipated that the same is true for the Schrodinger operator H — 
—A + V{x) with growing potentials like V{x) — |a;p™ treated in Theorem 12.21 



Theorem 13. II is obtained as a corollary from the following theorem of its inte- 
gral kernel version. Here one calculates explicitly the integral kernel K^'^\t, x,y) 
of [e~*^/^"e~*-^''/"e~*^/^"]" to estimate its difference from the integral kernel 
e~*^(a;,y) of e"*^. 



Theorem 3.2. (Azuma- Ichinose 2007 '2\) There exists a bounded operator R{t) 
and uniformly bounded operators {Q^"\t)}^^^-^ with integral kernels R{t,x,y) and 
Q^"'\t,x,y) being uniformy bounded continuous functions in (0,oo) x R x R such 
that 

K^'-\t,x,y)-e-'"{x,y) = [R{t,x,y) + Q^^\t,x,y)n'^]n'^- (3.2) 
they satisfy 

sup |i?(t, X, y)|, sup sup |(5'-"'(t, X, ?/)| ^ 0, t ^ ; sup |i?(i, x, ?/)| ^ 0, i ^ oo. 

a;,y n x,y x,y 

R{t,x,y) is explicitly given by 



R{t,x,y) = e--(x,y)^[t(i^;±^ + 



(e* + e *)xj/ — (x^ + y^) 



(e* — e 



t\2 



J_ / ^ 4xy-(e* + e~*)(x^+y^) ' 
16 V e* — e~* / 



(3.3) 



If t > 0, R(t,x,y) can become positive and negative. 
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Lemma 3.3. K'-''\t,x,y) 



1 + -^ 



1/2 



1 + 



X exp 



X exp-! 



2^/1 + 31. 



t n ( 
4n ~ 2t I 



1 



t2 , 

4n^ ^ 2r? 



1 " -a/1 + ^72+^ 



(3.4) 



Proof. Calculate the Gaussian integral 



n- 1 



f— 



R ./R 



n 

n[ 



dxi • • • (ia;„-i, 



where x — Xn, y — xq. We shall encounter with continued fraction to lead to the 
final expression (|3.4p of the lemma. 



To show Theorem 13.21 we simply calculate the difference K'-"'^ {t, x,y) — 
e~*^ {x,y), though it is not so simple. 

Here we mention what the operator with R{t, x, y) as its integral kernel is. By 
the Baker-Campbell-Hausdorff formula (e.g. [45], [10]), if A and B are bounded 
operators, we have 

jg-tB/2ng-tA/ng-tB/2njri _ ^-t{A+B) 

= exp ( - t{A + B)- n-^ — [2A + B, [A, B]] - Op(n-3)) 



+2 /.* 



24 7o 



-(t-s)(A+B) 



[2A + B, [A, B]]e-^(^+^)ds + Op{n-^), 



where Op(n''^) is an operator with norm of 0{n^^). In our case where A — 
— ^i9^, B = ix^, we can show R{t,x,y) is just the integral kernel of the oper- 
ator 

-'1 r. 

24 A 



-{t-s)H 



[2Ho + V,[Ho,V]]e-'''ds, 
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which does make sense, though Hq and V are unbounded operators. We have 
[2Ho + V, [Ho, V]] = -4Ho + 2V= -AH + W. 

4. Idea of Proof 

Put K{t) = e-^^/^e-^-^e-^^/^. Note that < K{t) < 1. Then we need to 
estimate the difference between K{t/n)" and e^*^ . The general technique of proof 
is: (i) to estabhsh an appropriate version of Chernoff's theorem ([S]): 

[(1 + r-i(l - K{r))-^ - (1 + H)-^] ^0, r i 

=> [Kit/nr - e-'"]^0, n^oo 

and/or (ii) to do telescoping: 

n 

fc=i 

to estimate each summand on the right. The former method (i) seems to be more 
efficient than the latter (ii). 

In fact, to prove Theorem 1, we use the former method, establishing the 
following norm version of Chernoff's theorem with error bounds. The case without 
error bounds was noted by Neidhardt-Zagrebnov [331 . 

Lemma 4.1. (Ichinose- Tamura [23 ) I. Let {i^(t)}t>o be a family of selfadjoint 
operators with < F{t) < 1, and H > a selfadjoint operator in a Hilbert space 
n. Define St := t-\l - F{t)). Then 

(a) For <a < 1, 11(1 + + tiO 
implies 

(b) For every fixed 6 > 0, 

|li^(i/n)"-e-*^|| =(5-2t-i+"e**0(n-"), n^oo, t > 0. 

Therefore for a — 1 this convergence is uniform on each compact interval 
[0, L] in the closed half line [0, oo). 

//. Moreover, in case H > rjl for some constant rj > 0, if for every e > there 
exists (5(e) > such that F{t) < 1 — (5(e) for all t > e, then 

||F(t/n)" - 6^*^11 = (l + 2/r/)2t-i+"0(n-"), n ^ oo, t > 0. 

Therefore for a = 1 this convergence is uniform on the whole closed half line 
[0,oo). 

Condition II is satisfied, e.g. for F(r) = e^'^^^'^e^'^^e^'^^/^. For the proof, 
we refer to [23] . 

For the proof of Theorems 2.2-2.5 we employ the latter method (ii), and 
further, for Theorems 2.2-2.4, make a crucial use of Agmon's kernel theorem: 
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Lemma 4.2. (Agmon's kernel theorem 

Let T : L'^{'R'^) L^(R'') be a bounded operator with ranges of T and its 
adjoint T* satisfying R[T], R[T*] C H"^{'R'^), m > d. If T is an integral operator 
with integral kernel T{x, y) being a bounded continuous function in R'^ x R*^ such 
that 

{Tf){x) = J Tix,y)f{y)dy, f ^ L\ 

then 

\T{x,y)\<C{\\TU + \\T*\U^\T\\^-^, 
where ||r||„j \\T\\c(l^^H"^) is the operator norm ofT as a bounded operator of 
L2(Rrf) into the Sobolev space W{R'^). 

Indeed, we estimate the £{L^ _ff'")-operator norm of the difference T = 

Jg-tV/2ngtA/ng-tV'/2nJn _ ^t{-A+V) 

5. Concluding Remarks 

We have so far considered the case where the operator sum H — A + B of two 
nonnegative selfadjoint operators A and B is selfadjoint. However, otherwise, the 
exponential product formula in norm does not in general hold for the form sum 
H = A + B oi two selfadjoint operators A > 0, B > 0, even if it is essentially 
selfadjoint on D[A] n D[B] (see 42i)- Nevertheless, there is some case where it 
holds: 

Theorem 5.1. (Ichinose-Neidhardt- Zagrebnov 2004 ^M) ^ ^ A+B be the 
form sum of A and B. If D[H°'] C DlA"] n ^[B"] for some ^ < a < I, and 
D[A^ C D[B^, then 

|j[g-tS/2ng-tA/„g^ti3/2n]„_g-tH|| ^ 0(„-(2"-l)), (5.1) 

||j^_tA/„g-tB/„j„_g-iff|j ^ 0(n-(2"-i)), (5.2) 

locally uniformly in [0, cxd). 

This error bound in (|5.1[) / (|5.2[) is also optimal. For this we refer to \4T . The 
condition for the domains of A and B is not symmetric. It is an open question 
whether one may improve it so as to become symmetric with respect to A and B. 

Finally, as we should like to mention, there is a very nice Feynman path 
integral formula which represents the Schrodinger semigroup, called the Feynman- 
Kac formula 

(e-*«/)(x) = (e-*(-^+^)/)(x) 
= / exp[- / V{B{s))ds]f{B{t))dti{B), 

J BeC{[0,oc)^R''),B{0)=x Jo 

where /i(-) is the Wiener measure on the path space C([0,oo) R'') (e.g. (Pj). 
We may use this formula to get whatever results, in fact, a lot of them. This is 



10 



Takashi Ichinose and Hideo Tamura 



a big advantage! But disadvantage is that it is only restricted to the Schrodinger 
operator or Laplacian. For instance, if we think of the semigroup for the relativis- 
tic Schrodinger operator H = \A^A + 77? + V{x), we have to estabhsh another 
Feynman-Kac formula (cf. p5]). 

Indeed, the Feynman-Kac formula is one of the realizations of Feynman path 
integral as a true integral on a path space. However, as Nelson |36j noted, the 
exponential/Trotter product formula also can give a meaning to the Feynman 
path integral as a time-sliced approximation by finite-dimensional integrals (cf. 
[15j). What it has advantage at is that we may apply it to the sum H = A + B oi 
any two selfadjoint operators B bounded from below. 
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